Current-Driven Switching in Magnetic Multilayer Nanopillars 
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We summarize our recent findings on how the current-driven magnetization switching in nanofabri- 
cated magnetic multilayers is affected by an applied magnetic field, changes of temperature, magnetic 
coupling between the ferromagnetic layers, variations in the multilayer structure, and the relative 
rotation of the layers' magnetizations. We show how these results can be interpreted with a model 
describing current-driven excitations as an effective current-dependent magnetic temperature. 



I. INTRODUCTION 



An exchange-based mechanism for current-driven 
switching of magnetization was predicted by Slon- 
czweski Q and Berger, and later observed experimen- 
tally. 0, 3 Recent research has been directed towards 
better understanding of the current-driven excitation 
mechanism, adequate description of the magnetic dynam- 
ics, and optimization of magnetic devices e.g. to decrease 
the switching current I s for possible ap plications in mag- 
netic memory, la El 0, H, 0, El illl Ha. fl4L Ua. Ua, H71 



We used the giant magnetoresistance effect (MR) to 
study current-driven switching in nanofabricated mag- 
netic F1/N/F2 trilayers (nanopillars). We studied how 
the current-driven switching is affected by variations of 
the magnetic field H, ambient temperature T p h, [2|| cou- 
pling between magnetic layers, 0i 113 electron spin- 
flipping in the spacer N or outside the trilayer, and mu- 
tual rotation of the two layers' magnetizations [3lJ • We 
interpret our experimental results for currents that are 
not too large in terms of the recently proposed effective 
temperature model, j^] 

Our samples were made with a multistep process 
described elsewhere. (24J Below, all thicknesses are in 
ran. The basic samples had structure Cu(80)/Fi(20- 
30)/N(10-15)/F 2 (2-6)/Cu(2-5)/Au(150). F was either 
Co or Py=Nig4Fei6. Although Co is commonly used 
in studies of current-driven switching, Py has the ad- 
vantage of small crystalline anisotropy and magnetostric- 
tion, allowing reproducible measurements at both 295 K 
and 4.2 K. The bottom Cu(80) layer was the extended 
lead, N, F2 and Cu(2-5) were patterned into an elongated 
shape with dimensions « (130 — 140) x (60 — 70) nm 2 , 
and Au(150) was the top lead. In all the samples, except 
for the AF-coupled ones (see Section. IlIIjL Fi was left ex- 
tended to minimize the effect of dipolar coupling on the 
current-driven switching. |Tlj Details of specific samples 
will be given with their data. We measured dV/dl with 
four-probes and lock-in detection, adding an ac current 
of amplitude 20-40 /iA at 8 kHz to the dc current /. Pos- 
itive current flows from Fi to F 2 . H is in the film plane 
and (except for the angular dependence studies) along 
the nanopillar easy axis. 
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FIG. 1: (a) Switching with current at 295 K in an uncoupled 
Py/Cu/Py trilayer. Solid line: H = 50 Oe, dashed line: H = 
— 500 Oe. Arrows mark the scan direction. It is the threshold 
current as defined in the text. Inset: MR dependence on H 
at / = 0. (b) Same as (a), at 4.2 K. (c),(d) Magnetization 
switching diagram, extracted from the current-switching at 
fixed H (solid symbols), and field-switching at fixed I (open 
symbols): (c) at 295 K, (d) at 4.2 K. Downward triangles: AP 
to P switching, upward: P to AP switching, circles: switching 
of the extended Py(20) layer. The reversible switching peaks 
are marked by coinciding upward and downward triangles. 
From Urazhdin et al. |25j| 



II. SWITCHING WITH VARIED H AND T, 



ph 



Figs. ^a,b) and insets show the variations of dV/dl 
with I and H for a Py(20)/Cu(f 0)/Py(6) sample with 
uncoupled Py layers at 295 K (la) and 4.2 K (lb), for 
H = 50 Oe (solid curves) and H = 500 Oe (dashed 
curves). The dependencies are similar to those seen in 
Co/Cu/Co samples. [H El At small H, the magneti- 
zation switches hysteretically to a higher resistance AP 
state at a large enough positive current I P ^ AP = I s > 0, 
and to a low resistance P state at lf p ^ p < 0. At larger 
H, the switching step turns into a, often much higher, 
nonhysteretic peak. The asymmetric /-dependence is a 
signature of the effect of the current on the magnetiza- 
tion, different from the Oersted field. Figs. IHc,d) show 
the switching diagrams at 295 K and 4.2 K, extracted 
from data such as those in Figs. Ea),(b), obtained both 
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FIG. 2: Differential resistance of a Py(20)/Cu(10)/Py(3.5) 
uncoupled trilayer at the marked values of H , at 295 K (left) 
and 4.2 K (right). Curves are offset for clarity. Dashed lines 
follow the reversible switching peak. 



by varying / at fixed H and H at fixed I. As expected, 
both the reduced magnetization and thermal activation 
result in smaller switching currents and fields H S (I = 0) 
at 295 K. The only other major difference between the 
295 K and 4.2 K data in Figs. |Ja,c) is the rounding of 
the 295 K hysteretic region at i < 0. In contrast, the 
hysteretic region at 4.2 K (Fig. QId)) is almost square. 

In the reversible switching regime, both the 295 K and 
4.2 K current-dependencies in Fig. ^ show a nearly lin- 
ear rise of dV/ dl above a threshold current It (labeled in 
Figs. m^jb)). In some samples, this rise displays struc- 
tures, clearly distinguishable from the switching peak due 
to their very weak dependence on H. Fig. shows data 
for a small Py(20)/Cu(10)/Py(3.5) sample, with dimen- 
sions estimated at 50 x 100 nm 2 . At 4.2 K, the linear 
rise is resolved into several peaks. They appear only to 
the left of the reversible switching peak, as the latter 
moves to higher I with increasing H (shown by a dashed 
line). At 295 K, the peaks are replaced by a smeared non- 
monotonic structure in dV j dl. Based on the data shown 
in Figs. ^ an d 121 and 9 other similar Py/Cu/Py samples 
studied, [25j we summarize the aspects of the threshold 

at 



P^AP 



behavior: i) I t is close to the small-H value of I. 
4.2 K. The correspondence I s « It is not universal, e.g. 
in magnetically coupled samples the switching may occur 
at I < (see Section ITTT)l . ii) Both I t and the sometimes 
observed structures in dV/dl at I > If are only weakly 
dependent on H, as compared to the reversible switch- 
ing peak. The 4.2 K data in Fig. El for H = 4 We 
give It ~ 1.0 mA, only 25% larger than I t ~ 0.8 mA 
at 250 Oe. iii) The height and inverse width of the re- 
versible switching peak in dV/ dl are correlated with the 
structure in dV/dl at / > I t : for example, in Fig. |2] (the 
295 K data), the reversible peak is tall and narrow when 
it is on the rising slope of the 'hump', but becomes wide 
and nearly disappears when on the trailing slope. 

We discuss first the nature of the threshold behavior 
starting at 7 t , and then the reversible switching peak at 
higher /. I t has been identified as the onset of large 
amplitude excitation of the patterned Py layer, as the 



magnetic energy provided by the current exceeds the lin- 
ear magnetic damping rate. |25| Starting at I t , there 
is a strong increase of excitation amplitude, and the 
highly excited magnetic state is manifested by increases 
in dV/dl, as illustrated in Figs. H HI I n an of om 
Co/Cu/Co samples, and most of our Py/Cu/Py samples, 
the excitations result in a nearly linear increase in dV/dl 
above It (see Fig. ^ and high H curves in Fig. Hfb)). 
In contrast, in Py/Cu/Py samples with small lateral di- 
mensions, at 4.2 K the excitations appear as a series of 
peaks in dV/dl above I t (Fig.[5J), corresponding to step 
increases in resistance V/I, with details that vary from 
sample to sample. We attribute these peaks to irregu- 
larities in sample shapes, resulting in a complicated de- 
pendence on I of both the magnetic damping rates and 
the distribution of excitations among different magnetic 
modes. In larger samples, the shape irregularities are less 
significant, giving a smoother dV/ dl. 

A different, quantum-mechanical origin has been pro- 
posed for the peaks, interpreted as the excitation thresh- 
old in point-contacts on magnetic multilayers. |j| This 
threshold would be due to the requirement for matching 
of the current-driven spin-accumulation with the low- 
est energies of the magnetic excitations, A// = hu>. In 
nanopillars, the quantum threshold value would double 
when H is increased from to « 1 kOe. Such a rapid in- 
crease is inconsistent with our data. In addition, even at 
4.2 K, thermal smearing with kT p h >> hcu w 10 jueV 
would completely smear out the quantum threshold. 
Thus, experiments rule out the quantum nature of the ex- 
citation threshold in nanopillars. We note that at large 
H > 2%M w 5 kOe, both the classical and quantum- 
mechanical thresholds depend linearly on H, and hui may 
become comparable to kT. In this limit, the quantum 
threshold may become important. 

The reversible switching peak in dV/dl is differ- 
ent from the threshold I t and the peaks we associate 
with magnetic excitations. Time resolved measurements 
showed that it is a consequence of telegraph noise with 
random distribution of dwell times in the P and AP 
states. 0; H3 Fig. GJa) shows the variations of aver- 
age dwell times tp(tap) in the P(AP) state with / for 
a Py(20)/Cu(10)/Py(6) sample at 295 K and 4.2 K. r P 
decreases as I increases, but tap increases. These vari- 
ations have a similar form at 295 K and 4.2 K, and are 
also similar to the variations in Co/Cu/Co. [l3| Fig.^b) 
shows that at a fixed I, Tp increases and tap decreases as 
| if | is increased (shown for 295 K). Fig.EJc) shows that, 
when both I and H are increased so as to hold Tp = tap , 
the average period of the telegraph noise decreases expo- 
nentially with similar slopes at 295 K and 4.2 K, down 
to the 1 MHz bandwidth limit of our setup. 

We now show how a peak in dV/dl at rp w tap can 
be derived from the variations of Tp, tap with /. For a 
fixed H, the average voltage across the sample is 



V(I) = I 



Raptap + RpTp 



Tp + TAP 



(1) 



3 




4.5 5.0 
I (mA) 



120 -115 
H (Oe) 




4 6 
l(mA) 



'(d) 








xaxr™ 




coxo 

1 


i i 







5 


10 



I (mA) 



FIG. 3: (a) The dependence of the P (downward tri- 
angles) and AP (upward triangles) dwell times on J, for 
a Py(20)/Cu(10)/Py(6) uncoupled sample. Open symbols: 
T ph = 4.2 K, H — -335 Oe, filled symbols: T ph = 295 K, 
H — —120 Oe. (b) Dependence of P (downward triangles) 
and AP (upward triangles) dwell times on H at / = 4.4 mA, 
T = 295 K. (c) Current dependence of the average telegraph 
noise period. H was adjusted approximately linearly with 
I, so that the average dwell times in AP and P states were 
equal. Solid circles: Py(20)/Cu(10)/Py(3.5) at 295 K, H =- 
93 to -121 Oe, open circles: same sample at 4.2 K, H =-300 to 
-315 Oe, solid squares: Py(20)/Cu(10)/Py(6) at 295 K, H = 
-113 to -123 Oe, open squares: same sample at 4.2 K, H =- 
300 to -450 Oe. (d) Circles: dV/dl vs. I at H = -0.3 kOe, 
4.2 K. Solid curve: a calculation, as described in the text, 
with To = 6.45 mA, and a + (3 = 27 mA -1 obtained from (a). 



with t p {I) k r exp[-a(7 - 7 Q )], t ap (I) rj r exp[/3(I 
To)], as in Fig. Eta)- We define 7 , r by T A p(h) 
Tp(Io) = To- From Eq. 



dV/dl « R P 
I (Rap — Rp) 



Rap — Rp 

l + exp[-(a + /?)(/- Jo)] 
exp[(a + /?)(J-Jo)](a + /3) 
(exp[(a + /3)(J-/ )]+l) 2, 



(2) 



The first two terms on the right are the resistance V/I, 
giving a step for the reversible transition from P to AP. 
For large lo (a + (3) , the last term has a maximum value 
Io(Rap — i?p)(cv + /3)/4 at I « Iq. This term gives rise to 
a peak in dV/ dl at I = Iq that can be much higher than 
Rap- Fig.Etd) shows a calculation (solid line) based on 
the data of Fig.Eta), and Eq. J3J), for I — 6.45 mA, and 
a + j3 = 27 mA -1 extracted from fig. Et a )- The calcula- 
tion is consistent with the dV/dl measurement (circles). 
We conclude that the reversible switching peak positions 
characterize the points (H, I) where rp = tap ■ The slope 
of the reversible switching line is a measure of the tele- 
graph noise variation with /, H. From Eq. El the height 



and inverse width of the reversible switching peak are 
given by a + (3, which is usually dominated by a, i.e. the 
dependence of the magnetic excitation rate in the P state 
on I. This is consistent with the correlation in Fig. [2Ja) 
between the reversible switching peak and the 'bump' 
associated with magnetic excitations: On the positive 
slope of the 'bump' (H = 150,170 Oe), the excitation 
level grows faster with I than on its the trailing slope. 
Consequently, a is larger on the rising slope of the bump, 
giving a tall and narrow reversible switching peak. The 
reversible peak nearly disappears on the trailing slope at 
H = 190, 210 Oe, where a is small. 

In contrast to I t , the switching peak does not repre- 
sent an onset of a physical process, it merely reflects the 
current-dependent telegraph noise statistics. The pres- 
ence of telegraph noise near the reversible switching line 
indicates that both AP and P states are unstable there. 
Strictly, the stability diagrams, Figs. dc,d) should be 
modified to include this unstable region. In most uncou- 
pled samples, the P-state becomes unstable at I ks It- 
The instability of the P-state at I > It is thus indirectly 
manifested in the rise of Rp. The AP-state is unstable 
at / both below and above the reversible switching peak. 
Since Tp is exponentially smaller than tap, the measure- 
ments of dV/ dl at I above the reversible switching peak 
give values very near Rap- 

As I and H are increased, rp, tap quickly approach the 
intrinsic nanosecond switching timescale (see Fig. Etc)), 
where telegraph noise is replaced by some sort of fast 
magnetic dynamics. Such a transition does not give any 
significant changes in the switching peaks in dV/dl (see 
e.g. Fig. [2| • Although the bandwidth limitation of our 
setup does not allow us to directly probe this regime, 
spectroscopic measurements show continuous spectrum 
from Hz up to « 1 MHz, characteristic of incoherent 
magnetic dynamics. |33| Microwave measurements show 
peaks at GHz frequencies on top of a broad background, 
evidence for some level of coherence. 281 



III. 



EFFECTS OF MAGNETIC COUPLING 



Fig. 0| summarizes the differences between uncou- 
pled (left), AF-coupled (middle), and F-coupled (right) 
Co(20)/Cu(d)/Co(2.5) samples at 295 K. For uncou- 
pled samples, we used d=10, and only patterned the 
top Co(2.5) and most of the Cu(10) layer to nanopil- 
lar size, leaving the bottom Co(20)-layer extended. F 
coupling was achieved by reducing the Cu thickness to 
Cu(2.6), near the third RKKY magnetoresistance (MR) 
minimum. 36] Dipolar AF coupling was achieved by also 
patterning about 10 nm of the Co(20) layer, with d = 6. 
Sample shape variations and interfacial roughness lead 
to variations in both AF and F coupling strengths. 

Fig. Ufa)) shows field-driven switching in uncoupled 
samples, similar to that in Fig.^for uncoupled Py-based 
samples. In AF-coupled samples (Fig.HJd)), the Co(2.5) 
layer is oriented AP to the extended Co(20) layer by the 
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FIG. 4: Results for uncoupled (a-c), AF-coupled (d-g), and 
F-coupled (h-j) Co/Cu/Co samples at 295 K. In the hysteretic 
plots, arrows show the scan direction. (a,d,h) dV/dl vs. H 
at 7 = 0. (b,e,i) dV/dl vs. I at the listed values of H. In (e) 
curves are offset for clarity. (c,f,g,j) Time-resolved measure- 
ments of R=V/I at the listed values of H, I. From Ref. |30l| 



dipolar field, so only the high resistance AP state is sta- 
ble at 77 = 0. In contrast to Fig. Ufa), both the P^AP 
and AP— >P transitions are now due to switching of the 
Co(2.5) layer. As 77 is scanned from a large negative to a 
large positive value, this layer switches three times: once 
to the AP state at negative 77, then together with the 
Co(20) layer at small positive 77 to stay in the AP state, 
and finally to a P state at large positive 77. The second 
transition sometimes produces a weak feature in dV/dl. 
In contrast, F-coupling causes the magnetizations to re- 
verse simultaneously at small H, giving only a small fea- 
ture in MR at 1=0 (Fig.^h)). In H-scans at large enough 
fixed I > 0, F-coupled samples gave 5% MR, similar to 
the 7 = values for uncoupled or AF-coupled samples. 

Fig.0Jb,e,i) compares the variations of dV/dl with 7. 
At small 77 = 20 — 50 Oe, applied to fix the magnetization 
of the bottom Co layer, the uncoupled sample (Fig.^b), 
solid curve) shows hysteretic switching, while the AF- 
coupled (Fig.^e), top curve) and F-coupled (Fig. Efi)) 
ones show reversible switching at V < and 7 > 0, re- 
spectively. At larger H, the switching in uncoupled sam- 
ples (Fig. Efb), dashed curve) becomes reversible, and 
in AF-coupled samples (Fig. E[e), lower two curves) it 
first becomes hysteretic and then reversible again. Time 
resolved measurements of resistance at 7, 77 near the re- 
versible switching peaks show telegraph noise switching 
between AP and P states (Fig.QJc,f,g,j)). At identical 7 
of opposite signs, the average telegraph noise periods in 
the AF-coupled sample are similar. 
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FIG. 5: Dependence of 1/I P ^ AP (upward triangles) and 
l/lf p ^ p (downward triangles) on AT? for uncoupled sam- 
ples. Open symbols: sample types 1 through 4, as labeled 
and explained in the text. Solid lines: best linear fits of data, 
excluding the angular dependence. The ordinate intercepts 
are zero within the uncertainty of the fits. 



IV. RELATION BETWEEN MR AND I s 

We studied the correlation between MR 
and I s in uncoupled samples with a structure 
Cu(80) /Py(30) /N(15) /Py(6) /Cu(2) /Au(150) . N=Cu 
in sample types 1 and 2. In sample types 3 
and 4, N was Cu(9.5)/Cu 94 Pt 6 (4)/Cu(1.5) and 
Cu(5.5)/Cu 94 Pt 6 (8)/Cu(1.5), respectively. The short 
295 K spin-diffusion length t s j » 6 nm (« 9 nm at 
4.2 K 35]) in Cu 96 Pt 6 decreases AR = Rap — Rp. 
We enhanced A77 in sample type 2 by replacing Cu(2) 
with a Cu(2)/Fe 5 oMn 50 (l)/Cu(2) sandwich. Fe 5 oMn 50 
is a strong spin-scatterer. j3f| Its placement between 
the trilayer and the top lead reduces the negative effect 
of the spin accumulation outside the trilayer on the 
MR. We also examined the variation of AT? and the 
switching currents 7 S in samples of type 1 while rotating 
the magnetization of the Py(30) with a small 77 in the 
film plane (similar results were shown in (2(J). 

The results, averaged over at least 7 samples of each 
type, are summarized in Fig. where data for angular 
dependence are also shown. Both lf p and lf p follow 
an approximately inversely linear dependence on AR, re- 
gardless of the method by which A77 was varied. 



V. EFFECTIVE TEMPERATURE MODEL 

In their pioneering work, Slonczweski [l| and 
Berger, predicted that exchange interaction leads 
to current-driven magnetic excitation by spin-polarized 
current. Berger considered generation of magnons by 
spin-flipping of electrons, driven by the spin accumu- 
lation. Slonczewski considered the electron spin com- 
ponent transverse to the magnetization, which is ab- 
sorbed by the ferromagnet due to a combination of spin- 
dependent reflection at the interfaces and averaging of 
the spin precession phases in the ferromagnet. The re- 



suiting torque drives the magnetic dynamics. This model 
postulates conservation of the total magnetic moment, 
and thus captures only the coherent (uniform) magnetic 
dynamics. It can not consistently treat the excitation 
of finite-wavelength spin- waves, whose generation is gen- 
erally not associated with transfer of angular momen- 
tum transverse to the magnetization (a similar argument 
is used to prove that only uniform precession is excited 
in the transverse FMR experiments). Similarly, Berger 
made the approximation that only the uniform preces- 
sion is generated by electron spin-flipping, leading again 
to coherent magnetic dynamics. 

Since the characteristic magnon energies are signifi- 
cantly lower than the typical conduction electron ener- 
gies (see the above discussion of the lack of quantum 
threshold behavior), we assume that a large number of 
magnetic modes are excited by the current. The popula- 
tions rii of the modes with energies Ei can be then ap- 
proximately described by a single parameter, an effective 
temperature T m . so that rii « kBT m /Ei for the degen- 
erate modes. [2J, |25|, |32| This approximation fails when 
T m approaches the Curie temperature of the ferromag- 
nct. We emphasize that this is just a single-parameter 
approximation for a generally much more complex, non- 
thermalized excitation distribution. Because of the large 
magnon populations, the spontaneous (independent of 
rii) magnon emission can be neglected. In a ballistic 
transport approximation [3^ | 



kT n 



kT, 



pi) 



\ + 2peVB/ 1 ' 



(3) 



where B is a constant characterizing the strength of the 
exchange interaction, 7 is related to the Gilbert damping 
parameter in the classical Landau-Lifshitz equation, V 
is the voltage across the trilayer, and p is the current 
polarization, created at the location of F2 by Fi, if F2 is 
removed, p > 0(< 0) in the P(AP) state for Py/Cu/Py or 
Co/Cu/Co trilayers, and p = if Fi is removed. Eq.|3]is 
similar to the expression obtained by Berger in a diffusive 
transport approximation (Eq. (9) in Ref. @). 

Eq. J2Jl diverges at V — » -j/(2peB). We identify this 
divergence with the threshold It for the large amplitude 
of excitations. Eq. © is not applicable above the thresh- 
old, because the conduction electrons are scattered with 
a large spin-flip probability (« 0.5), nearly independent 
of T m . For this regime, we have proposed an empirical 
relation Hal 



Tin. 



ph 



K{I-I t ) for / > I t 



(4) 



where if is a constant determined by the magnetic re- 
laxation rate. In the thermal activation model, the dwell 
times tp^ap are determined by 



I 

Tp,AP = Q eX P 



Up 



AP 



kT, 



P,AP 



(5) 



where Q « 10 7 s 1 [^ij is the effective attempt rate, 
Up t AP is the potential barrier for switching from the P or 




FIG. 6: Schematics of current-driven switching, as explained 
in the text. Dashed lines indicate T m . 



AP state. T m is approximated by Eq. J3J| for the low-H 
hysteretic switching, or Eq. |(2J for the telegraph noise. 
The areas of the P(AP) state stability are determined 
by Tp(AP) > f s. Up^ap depend on I only through the 
variation of the magnetization with T m . 

Fig. HJa) schematically describes hysteretic current- 
driven switching in uncoupled samples at small H. From 
Eq. 0, T£ drastically increases when I approaches It, 
giving rise to a thermally activated transition into the 
AP state at JcT£ » ,_,V p „ s ~ for the data ac- 



ln(t cxp n) ~ 1 

quisition time t exp = I sec/point. In the AP state, 
Eq. © predicts that the magnetic system weakly cools 
to T^ p < T p h, so that at H < H s it becomes trapped 
in this state. Fig. HJb) is a schematic for H > H s , giv- 
ing kTph > ^ff-. The AP to P switching is now ther- 
mally activated, resulting in telegraph noise at / > I s , 
H > H s . Fig.|SJc) is the schematic for AF-coupled sam- 
ples at small H . AF coupling increases Uap and reduces 
Up, so that P to AP switching is now thermally acti- 
vated. At large enough I < 0, reverse switching also 
becomes thermally activated, giving telegraph noise. As 
H grows, the schematic for AF-coupling becomes 6(a), 
giving hysteretic switching, and then 6(b), giving tele- 
graph noise for large enough positive /. Finally, 6(b) is 
the schematic for F-coupling at all H, so that AP to P 
switching is always thermally activated and large enough 
I > gives telegraph noise. 

The difference between the 295 K and 4.2 K switching 
diagrams (Fig. is qualitatively well described by the 
effective temperature model. |32j The rounding of the 
295 K diagram at / < is due to enhancement of ther- 
mal fluctuations of magnetization even at small I (see 
Eq. |2l . Increasing H decreases Uap , so that even these 
weak excitations can activate the AP— >P transition. In 
contrast, the almost square 4.2 K diagram reflects the 
small excitation amplitude at I < I t , beyond which the 
excitation rapidly grows, giving the switching. 

We use Eqs. (@J and |J5J to describe the dependen- 
cies of Tp on /, H, and T p h in Fig. [21 Eq. © should 
be used for tap(I, H, T p h). These equations capture the 
tendencies qualitatively well. From the slope of Tp(I), 
the increase rate of in Py(20)/Cu(10)/Py(6) samples 
is estimated at « 400 K/mA. [29] The thermal activation 
model breaks down when kpT m becomes similar to Up. 
This limitation of the model should be understood in the 
context of the actual magnetic dynamics. T m ~ U in 
Eq. [3[ gives a switching rate similar to the intrinsic mag- 
netic dynamic rates. Telegraph noise switching is then 
replaced with fast fluctuations of the magnetic moment. 
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The effective temperature concept (denned for a given 
magnetic orientation) becomes irrelevant in this regime. 

We note that, although the 295 K and 4.2 K current- 
dependencies in Fig. EUa) are given for different values 
of H, the effect of this difference on rp should be small, 
as can be seen from Fig. Etb). Since Up is large, its 
relative variation with H is small, compared to much 
larger relative variation of Uap- Thus, Fig. EJa) can be 
used to determine Tp(I), while tap is adjusted with H. 

Finally, to describe the inverse relationship in Fig. [S] 
we use the simplest plausible model, in which AR is pro- 
portional to the current polarization p at the location of 
F 2 due to Fi. This crude approach is self-consistent: the 
MR is correctly predicted to disappear if Fi is absent. 
The current-driven switching occurs at 16fc_BT m rs U. 
From Eq. [21 follows 1/ I s oc p, giving the dependence in 
Fig. In this context, the angular dependence can be 
understood similarly in terms of the projection of p onto 
the direction of the magnetization of F 2 . 

VI. SUMMARY 

We summarize the following important experimen- 
tal observations for current effects in magnetic trilayer 



nanopillars: i) a square switching diagram at 4.2 K, and 
rounded at 295 K; ii) an onset current I t (closely related 
to the hysteretic switching current I s ) for a linear rise of 
dV/dl in larger samples or a series of peaks in smaller 
ones; iii) reversible switching, characterized by telegraph 
noise with rate both increasing exponentially with I and 
shifting with temperature. The reversible switching peak 
in dV/dl occurs when the dwell times in the P and AP 
states are approximately equal, iv) Reversible switching 
at small H, if magnetic layers are F- or AF-coupled. v) 
Inverse linear relation between MR and I s . We explain 
the observed behaviors in terms of thermal activation 
over a magnetic barrier, with a current driven effective 
magnetic temperature. 
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